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The curvature field is measured from tracer particle trajectories in a two-dimensional fluid flow 
that exhibits spatiotemporal chaos, and is used to extract the hyperbolic and elliptic points of 
the flow. These special points are pinned to the forcing when the driving is weak, but wander 
over the domain and interact in pairs at stronger driving, changing the local topology of the flow. 
Their behavior reveals a two-stage transition to spatiotemporal chaos: a gradual loss of spatial and 
temporal order followed by an abrupt onset of topological changes. 



PACS numbers: 47.52. +j, 47.20.Ky, 05.45.-a 

When a system governed by nonlinear equations of mo- 
tion is driven out of equilibrium, a variety of complex be- 
haviors can result, ranging from chaos in low-dimensional 
systems [T] to the seemingly random dynamics of turbu- 
lent flows [2] . If the number of active degrees of freedom 
(roughly corresponding to the number of equations re- 
quired to characterize the dynamics) is small, the math- 
ematics of dynamical systems has proved to be a power- 
ful tool. When it is quasi-infinite and the dynamics are 
turbulent, statistical approaches based on assumed scale 
invariance have been fruitful. For spatially-extended sys- 
tems below the turbulence transition, however, a regime 
with a large but finite number of active degrees of free- 
dom exists that is disordered both in space and in time 
[U m |SJ ini 13 [SI [U . This regime of spatiotemporal chaos 
(STC) remains poorly characterized and understood; in- 
deed, there is no generally agreed-upon quantitative in- 
dicator of STC. 

In this Letter, we study the dynamics of a simple 
fluid flow that exhibits STC by considering its under- 
lying topology. We describe a method for locating the 
time-dependent topologically special points of the flow, 
and show that their dynamics describe the flow pattern 
as a whole. These special points undergo pairwise inter- 
actions, changing the flow topology when they are cre- 
ated or annihilated. Surprisingly, we flnd that the rate 
of creation or annihilation shows a discrete onset, while 
other measures of spatial and temporal disorder increase 
smoothly with the rate of energy input. This form of 
STC is characterized not only by disorder but also by 
constantly changing flow topology. 

In a driven body of fluid, there can exist instantaneous 
stagnation points where the velocity vanishes, relative to 
some observer. These special points carry the bulk of the 
information contained in the flow fleld: if the locations 
of all of these topologically special points and their local 
flow properties are known, most of the full flow fleld can 
be determined [TU] . They are also an essential component 
of chaotic mixing These points are distinct from 

the topological defects previously considered in studies 
of STC [m |T3J [13] , as they are present even when the 



flow is not chaotic. In a two-dimensional (2D) flow field, 
the special points come in two types. When embedded 
in a region of the flow that is dominated by vorticity, 
they are elliptic; in a strain-dominated region, they are 
hyperbolic {i.e. saddle-like). These special points, how- 
ever, have proved to be very difficult to identify, particu- 
larly in experimental flows. In this Letter, we show that 
by considering the curvature of Lagrangian trajectories, 
that is, the trajectories of individual moving fluid ele- 
ments, we can flnd the elliptic and hyperbolic points in 
an automated way, even when they move. Once located, 
the trajectories and statistics of the special points give 
insight into the transition to and dynamics of STC. 

We generate a quasi-2D flow using magnetohydrody- 
namic forcing in a thin layer of conducting fluid, as 
described previously |15L I16j . A 6 mm layer of water 
containing 8% by weight of CUSO4 was placed above 
a square lattice of permanent magnets with alternating 
orientation. When a current is driven across the cell, 
Lorentz forces set the fluid into motion. The dimension- 
less strength of the forcing is measured by the Reynolds 
number Re = UL/v, where U is the root-mean-square 
velocity, L = 2 cm is the mean magnet spacing, and 
1/ is the kinematic viscosity. At low Re, the flow is a 
square array of vortices of alternating sign, as shown in 
Fig. 1. As Re grows, however, the flow deviates from 
the forced lattice and becomes spatiotemporally chaotic. 
To measure the flow, we follow the simultaneous trajecto- 
ries of thousands of neutrally-buoyant 116 /im fluorescent 
polystyrene tracer particles, using algorithms similar to 
those described by Ouellette et al. [T7]. The particles are 
imaged at a rate of 12 Hz, and their positions are deter- 
mined to a precision of 25 fim (0.1 pixels). The velocities 
and accelerations of the particles are then computed by 
fitting polynomials to short segments of the trajectories 
[T5] . Statistics are collected in a 10 cm x 10 cm window 
in the center of the flow, so that boundary effects are 
excluded. 

To find the elliptic and hyperbolic points, we first com- 
pute the curvature along the trajectories of the tracer 
particles. Curvature, a geometrical quantity containing. 
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FIG. 1: (color online) Velocity, vorticity, and curvature fields for Re = 32 (a, d), 93 (b, e), and 245 (c, f). In (a-c), velocity 
vectors are shown as arrows, undersampled by a factor of 8 for clarity. The vorticity is shown by color: red (or gray) corresponds 
to large negative vorticity (clockwise rotation), and blue (or black) to large positive vorticity (counterclockwise rotation). As 
Re increases, the flow becomes more disordered, (d-f) show the logarithm of the curvature; red corresponds to large and blue 
to small values. Low values of curvature typically form lines, while the highest values appear as points. 



in principle, no dynamical information, completely spec- 
ifies a curve in 2D space. Because of the nature of our 
measurement technique, however, the trajectories are pa- 
rameterized by time. In this case, the Frenet formulas 
show that the curvature is given by k = a„/u^, where a,i 
is the acceleration normal to the direction of motion and 
u is the velocity of the particle. The single-point statis- 
tics of curvature have previously been studied in turbu- 
lent flow [THJ , but were found to be explainable with 
a simple model that should also apply to non-turbulent 
flows. In our flow, we measure curvature probability den- 
sity functions (PDFs) that are consistent with the model 
proposed by Xu et al. [20] . 

Instead of focusing on such single-point curvature 
statistics, we consider curvature fields, analogous to ve- 
locity or vorticity fields. Sample curvature fields for the 
steady and STC regimes are shown in Fig. 1, and strik- 
ing structure is evident. As shown previously in studies 
of turbulence, the distribution of curvature is exception- 
ally wide [m [20]. What was not observed in previous 
studies, however, is the tendency of low values of curva- 
ture to be spatially organized into lines, while high values 
exist as solitary points. Comparing the vorticity and cur- 
vature fields in Fig. 1, we see that these high-curvature 



points correspond to the hyperbolic and elliptic points 
of the flow. This observation has a clear physical inter- 
pretation: near both hyperbolic and elliptic points, the 
direction of fluid particle trajectories changes over very 
short length scales, corresponding to intense curvature. 
By locating the local maxima (with values larger than 
the mean) of the curvature field, we can therefore find the 
topologically special points of the flow. To classify them 
as elliptic or hyperbolic, we make use of the Okubo- Weiss 
parameter A = (w^ — s'^)/4, where is the enstrophy 
and is the square of the strain rate [HI [23] . If a curva- 
ture maximum lies in a region with A > 0, where rotation 
dominates the flow, we classify it as elliptic; if A < 0, the 
local flow is dominated by strain and we classify the point 
as hyperbolic. 

Once we have located the hyperbolic and elliptic points 
at each instant in time, we can feed their positions into 
the same tracking algorithms we use to construct the 
tracer-particle trajectories, and thereby study their dy- 
namics. At low Re, where the underlying flow is a vortex 
array locked to the magnetic forcing, the hyperbolic and 
elliptic points lie on a square lattice, as shown in Fig. 2, 
with the elliptic points in the vortex centers and the hy- 
perbolic points at the vortex corners. As Re is increased. 
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FIG. 2; (color online) Trajectories of the topologically special 
points for the three Reynolds numbers in Fig. 1. Hyperbolic 
points are plotted in red, while elliptic points are in blue. 
At Re = 32 (a), the special points are tightly bound to the 
forced vortex lattice. At Re = 93 (b), they remain bound but 
make larger excursions. At Re — 245 (c), where the flow is 
spatiotemporally chaotic, the special points wander over the 
domain. 



the special points wander progressively farther from their 
lattice sites. We observe that the elliptic points move in 
quasi-circular orbits, while the hyperbolic points move 
primarily along their stable manifolds. Finally, when Re 
is high enough, the special points break free from the 
lattice and wander freely; in a sense, the lattice melts. 
We note that the dynamics of the hyperbolic and ellip- 
tic points are representative of the dynamics of the flow 
pattern, rather than of the tracer particles themselves; 
the individual tracer particles may wander between the 
vortex cells even in the regime where the special points 
remain pinned. 

To elucidate the transition between the ordered and 
disordered states of the topologically special points, we 
show in Fig. 3 the radial distribution function g{r) of 
the special point positions for three Reynolds numbers. 
g{r) is defined to be the mean number density a distance 
r from a fixed position, normalized by the bulk number 
density. In an ordered state, g{r) is expected to show a 
series of peaks corresponding to lattice sites; in a disor- 
dered state, however, g{r) should be unity. At Re — 32, 
when the special points are tightly bound to their lattice 
sites, g{r) is found to be peaked at many of the loca- 
tions expected for a 2D square lattice. As Re increases, 
the peaks broaden and g{r) gradually loses structure. By 
Re = 245, where the special points move freely and the 




r/L 



FIG. 3: (color online) The radial distribution function g{r) 
of the topologically special points, shown for Re — 32 (□), 93 
(O), and 245 (A). The separation r is scaled by the mean 
magnet spacing L. Dashed vertical lines show the positions 
of some of the peaks expected for a 2D square lattice (other 
peaks are not observed due to the finite resolution). As Re 
increases, the spatial order vanishes and g{r) becomes liquid- 
like. 




FIG. 4: (color online) The maximum value of g{r) (□, left 
axis), showing a gradual transition to spatial disorder with 
increasing Re, and the annihilation rate (O, right axis) of 
hyperbolic-elliptic pairs, which shows a well-defined thresh- 
old. The solid lines are drawn to guide the eye. 



spatial order has vanished, g(r) is unity and the special 
points form a liquid- like state. To quantify the loss of 
order, we plot the maximum value of g(r) (correspond- 
ing to the height of its first peak) as a function of Re in 
Fig. 4. 

Once the topologically special points can move appre- 
ciably around their lattice sites, they undergo pairwise in- 
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FIG. 5: (color online) PDFs of the number A*' of topologically 
special points present in the measurement area at each instant 
in time, normalized by the number of lattice sites, for four 
Reynolds numbers. As Re increases, the mean number of 
special points decreases, while the width of the distribution 
grows. 



teractions that change the topology of the flow. They can 
be annihilated in vortex mergers, and new hyperbolic- 
elliptic pairs are created when vortices split. By record- 
ing the number of such events from our special point 
trajectories, we can measure the rates for each of these 
processes. The annihilation rate is shown as a function 
of Re in Fig. 4; we find that the creation and annihi- 
lation rates are equal to within experimental accuracy. 
As expected, these rates grow substantially as the driv- 
ing increases. Surprisingly, however, these rates have a 
clear onset at a critical Re, unlike the gradual decline 
of g{r). We have also measured other indicators of dis- 
order, including the average of the local speed (du/dt) 
and the mean of the largest Lyapunov exponent (A). As 
with g{r), however, these measurements show a continu- 
ous change rather than a sharp onset. The lack of a clear 
threshold in these measurements may be attributed to 
the averaging of g{r), {du/dt), and (A) over space, time, 
or both. In contrast, the creation and annihilation of hy- 
perbolic and elliptic point pairs are local in both space 
and time. 

The number of special points in the measurement vol- 
ume changes in tandem with the rise of the pair creation 
and annihilation rates. In Fig. 5, we show the PDFs of 
the number of special points at each instant in time for 
several Reynolds numbers. As Re increases, the pattern 
coarsens and the mean number drops, consistent with the 
well-known inverse energy cascade in 2D turbulence [53] . 
At the same time, the width of the distribution grows, 
which signals the increased activity of the special points 



and the corresponding increase in the frequency of topo- 
logical changes. 

In summary, we have developed a method to locate 
the moving hyperbolic and elliptic points experimentally 
by measuring the curvature of particle trajectories, and 
have used them to characterize the dynamics of a 2D flow. 
As Re is increased, these points gradually unbind from 
their preferred locations (determined by the forcing) , de- 
forming the entire flow pattern. When they are created 
or annihilated in pairs, starting at Re ^ 75, they change 
the topology of the pattern. The behavior of these points 
shows that the transition to STC involves two successive 
stages: a gradual loss of spatial and temporal order fol- 
lowed by a surprising abrupt onset of topological changes. 
We suggest that new theories of STC may be developed 
using these topological ideas. 

We thank D. Egolf, P. Love, H. Riecke, and G. Voth 
for helpful discussions. This research was supported by 
the U.S. National Science Foundation under Grant DMR- 
0405187. 
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